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ABSTRACT 

Intergalactic absorbers along lines of sight to distant quasars are a powerful diagnostic for the 
evolution and content of the intergalactic medium (IGM). In this study, we use the FUSE satellite to 
search 129 known Lya absorption systems at z < 0.15 toward 31 AGN for corresponding absorption 
from higher Lyman lines and the important metal ions VI and C III. We detect VI in 40 systems 
over a smaller range of column density (logVovi = 13.0-14.35) than seen in H I (logVni = 13.0- 
16.0). The co-existence of VI and H I suggests a multiphase IGM, with warm photoionized and 
hot ionized components. With improved VI detection statistics, we find a steep distribution in 
VI column density, dAfovi/dNovi oc VQy T 2±01 , suggesting that numerous weak O VI absorbers 
contain baryonic mass comparable to the rare strong absorbers. Down to 30 mA equivalent width 
(O VI A1032) we find an absorber frequency dNowi/dz m 17 ±3. The total cosmological mass fraction 
in this hot gas is at least ^ W him = (0.0022 ± 0.0003)[/i 7 o(2'o/0.12 , o)(/ovi/0.2)]- 1 where we have 
scaled to fiducial values of oxygen metallicity, O VI ionization fraction, and Hubble constant. Gas in 
the WHIM at 10 5-6 K contributes at least 4.8 ± 0.9% of the total baryonic mass at z < 0.15. We 
then combine empirical scaling relations for the observed "multiphase ratio" , Vhi/Vovi oc Vm ' 9 ^' 1 , 
and for hydrogen overdensity in cosmological simulations, A/hi oc <5^ 7 , with the H I photoionization 
correction to derive the mean oxygen metallicity, Zq w (O.O9Z0)(/ovi/O.2) _1 in the low-z multiphase 
gas. Given the spread in the empirical relations and in /ovi, the baryon content in the O VI WHIM 
could be as large as 10%. Our survey is based on a large improvement in the number of O VI absorbers 
(40 vs. 10) and total redshift pathlength (Az ks 2.2 vs. Az « 0.5) compared to earlier surveys. 
Subject headings: cosmological parameters — cosmology: observations — intergalactic medium — 
quasars: absorption lines 



1. INTRODUCTION 

One of the great, unanticipated legacies of the Far Ul- 
traviolet Spectroscopic Explorer (FUSE) mission is surely 
the detection of O VI absorption lines along extragalac- 
tic sight lines from the warm-hot intergalactic medium 
(WHIM) at T = 10 5 " 7 K, as predicted by cosmolog i- 
cal simulations ijCen fc Ostr ikcr 19991 iDave et al] l2001L 
Hot gas in the intergalactic medium (IGM) is produced 
by shocks generated by gravitational instability during 
the formation of large-scale structure, and its detection 
in highly ionized oxygen is an indicator of widespread 
metal transport into the IGM through feedback from 
galaxy formation. The O VI absorption probes gas at 
10 5 " 6 K, somewhat cooler than the bulk of the WHIM, 
which has been shock-heated to temperatures of several 
xlO 6 K. This latter, very hot gas is only detectable 
through weak X-ray absorption lines from higher ions 
(e.g., O VI I O VIII, Ne V I II, N VI, N VII) as described 
recently bv lNicastro et al.l ()2005[) . 

The search for the WHIM gas has now begun in 
earnest, using sensitive UV resonance absorption lines, 
primarily the O VI doublet (1031.926, 1037.617 A) which 
is ~100 times more sensitive than the X-ray transi- 
tions of O VII (21.602 A) or O VIII (18.97 A). Thus, 
the initial WHIM detections of O VI were made by 
FUSE at z < 0.15 and by the Hubble Space Telescope 
(HST) at z > 0.12, but over a modest total redshift 



pathlength Az » 0.5 (jTripp. Savage, fc Jenkmsl 1200(1 
iSavaee et aill2002j) . For the portion of WHIM contain- 
ing O VI, the estimated baryon fractions were ~5% 
for O/H metallicity equal to 10% of th e former solar 
value, 7.41 x 10~ 4 ijGrevesse et al.l[l996|) . Because the 
solar oxygen abundance has recently been re-measure d 
at (O/H) = 4.90 x 10" 4 (lAllende Prieto et al.ll200l . 
we have re-scaled our baryon estimates to this lower 
value. These estimates have large uncertainties, includ- 
ing untested assumptions of uniform metallicity, O VI 
ionization equilibrium, and multiphase IGM structure. 

The initial O VI searches r eported only a handful of 
absorb ers - four systems bv iTripp. Savage. &: Jenkinsl 
(2000) and six bv lSavage et all l|2002|) . At low redshift. 
down to 50 mA equivalent width in O VI A1032, the 
number of O VI absorbers per unit redshift is 10-15% 
of that found in the H I Lya surveys: dNovi/dz w 
14 + 9 JSavage et all l200l versu s dM m /dz « 112 ± 
9 (jPenton. Stocke. fc ShuTH I2004J) . In this study, we 
dramatically increase the total surveyed pathlength to 
Az > 2 and the total number of O VI absorbers to 
40. We also study the relationship between absorbers 
from the WHIM and the warm neutral medium (WNM; 
10 3 - 5-4 - 5 K). With the increased number of absorbers, 
we can begin to look at WHIM absorption in a statisti- 
cal manner, particularly the distribution in O VI column 
density. 



1 also at JILA, University of Colorado and National Institute of 
Standards and Technology 



2. OBSERVATIONS AND DATA ANALYSIS 
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Fig. 1. — Comparison of HI and OVI column densities in 124 known Lya absorbers. The left panel shows a mild correlation, if any, 
between Nm and Novi- While TVjji varies over a factor of nearly 1000, iVovi varies by a factor of only 20. The right panel shows the 
"multiphase ratio", iVjn/ iVoyij with a dispersion of less than one dex over nearly three orders of magnitude in Nm, suggesting that 
photoionized H I and warm-hot gas occupy different phases of the IGM. The multiphase ratio is well fitted by a power law, -/Vjji/ A^ovi 
= 2.5 X [Njji/10 14 cm -2 ] ' 1 (dashed line). TVhi can be related to the overdensity, Su = p/{p), as discussed in § 4 and shown along the 
top axis. 



We began our survey with published lists of 
intervening Lya absorption systems toward low-2 
AGN obtained from GHRS and STIS surveys by 
IPenton. Shull. fc Storied Eml : IPenton . Stor.ke. fc Shulll 
(2004). Other sight lines were covered by literature 
sources or measured at the University of Colorado from 
STIS/E140M spectra as discussed in the upcoming Dan- 
forth, Shull, & Rosenberg (2005, hereafter Paper II). We 
disregarded any weak Lya absorbers (Wa < 80 mA, 
log TVhi < 13.2) and searched the FUSE data for higher- 
order Lyman lines plus O VI and C III counterparts. We 
have been conservative in our identification of O VI sys- 
tems, requiring unambiguous (> 4<r) features, and exam- 
ining multiple channels of FUSE data and both lines of 
the O VI doublet when possible. In all, we analyzed 171 
absorbers in 31 FUSE sight lines. We measured all avail- 
able Lyman series lines to determine accurate Nm and 
doppler b values for H I v ia curve-of-growth concordance 
curves (|Shull et alJl2000l) . 

Of the Lya absorbers, 129 were at z < 0.15, where 
O VI absorption could potentially be observed. We de- 
tected O VI at 4cr or greater level in one or both lines 
of the doublet for 40 absorbers, and we obtained 4er or 
greater upper limits in 84 other cases. The remaining 
five absorbers fell on top of airglow or strong lines from 
the interstellar medium (ISM), were blended, or were in 
some other way inaccessible. 

We determined metal-ion column densities via Voigt 
profile fits and /or the apparent column method 
ijSavaee k, Sembachlll991|) . We assumed that any satu- 
ration in these lines was mild, and that profile fits accu- 
rately determine the column density. In cases where we 
detected both O VI lines, we assigned a weighted mean 
of the column densities to Nqvi- We describe the full de- 
tails of absorber selection and data analysis in Paper II. 



3. THE MULTIPHASE IGM 

The O VI lines are ideal tracers (at io 5 5±0 - 3 K) for 
the warm- hot ionized medium (WHIM; 10 5 - 10 7 K), 
while H I traces the photoionized warm neutral medium 
(WNM; 10 3 - 5 " 4 - 5 K). The left panel of Figure 1 shows a 
weak correlation, if any, between iVovi and Nja- While 
./Vhi varies over nearly three orders of magnitude, iVovi is 
detected between 10 13 cm -2 and a few times 10 14 cm~ 2 
(a factor of ~20) . To gauge the relative amounts of warm 
photoionized gas (H I) and hot collisionally ionized gas 
(O VI), we use the "multipha se ratio", N m/NgvT. This 
ratio was defined previously ijShullI 12003) as a means of 
assessing the range in contributions from photoionized 
gas (H I) and collisionally ionized gas (O VI) , which often 
appear to be associated kinematic ally ijSembach et alJ 
l2003tlCollins. Shull fc Girouxll2004]l . 

As shown in the right panel of Figure 1, the Nju/Nqvi 
ratio exhibits a strong correlation with Nm, with a typ- 
ical range Nm/Novi ~ 0.5-50 and a dispersion of less 
than one (dex). Higher H I column absorbers typically 
display a higher multiphase ratio, while weak H I ab- 
sorbers have iVovi w (0-1 — 2) ./Vhi. We are well aware of 
the fact that plots such as Figure 1 are subject to the ef- 
fects of correlated errors in Nm- However, this ratio has 
physical utility in deriving a mean (O/H) metallicity in 
multiphase gas (see § 4) , when combined with the empir- 
ical correlation between TVhi and hy drogen overdensit y, 
Sh, seen in cosmological simulations ( D ave et al.H l999). 

The large range in TVhi and the correlation in multi- 
phase ratio provide evidence that the IGM has at least 
two phases (WHIM and WNM). If the H I and O VI 
materials were well mixed, we would expect that the 
multiphase ratio would be constant with -/Vhi. We sug- 
gest instead that the WHIM occupies a shell around a 
warm neutral core. The outer WHIM shell is heated by 



a combination of external ionizing photons (from AGN) , 
shocks from infalling clouds, and possible shocks from 
cluster outflows such as superwinds and SNR feedback. 
The narrow range of iVovi compared to Nm implies that 
the WHIM shell may have a characteristic column den- 
sity of 10 13 — 10 14 cm -2 , while the neutral core can be 
arbitrarily large. 

In this scenario, the very large values of the multi- 
phase ratio may arise in "quiescent gas" , possibly lo- 
cated in high-column density H I gas in halos. The 
absence of associated O VI suggests the lack of shocks 
at velocities greater than about 150 km s _1 . Many 
strong H I (Lya) absorbers have been seen in proximity 
(with i n 200/i^n 1 kpc) of bright galaxies jLanzet 
fl995t IChen etalJ 11993 IPenton. Stocke. fc Skill 1200? 
In our current O VI sample, we found three absorbers 
with extremely high multiphase ratios: very high neu- 
tral hydrogen columns with no detectable WHIM. In the 
first case, the absorber at cz — 1586 km s _1 toward 
3C273 shows \ogN m = 15.67 and logiVbvi < 13.17 and 
is 71/iyQ 1 kpc on t he sky from a dwarf {Mb — —13.9) post- 
starburst galaxy ijStocke et al. 2004) . The second case is 
the absorber at cz = 23,688 km s _1 toward PHL1811 
with \ogN m = 15.94 and logTVovi < 13.06. This sys- 
tem is beyond the L* survey depth, so we cannot com- 
ment on surrounding galaxies. Finally, the Lyman-limit 
system at cz — 24,215 km s _1 toward PHL1811 shows 
logiVHi = 18.11 and log-ZVovi < 13.06. This absorber 
lies 23" {Mh^ kpc) fr om an L* g a laxy a t z = 0.0808 
I) Jenkins et al.l 12003^ . iTripp et all $200!$ note a sub- 
DLA system at z abs = 0.00632 toward PG 1216+069 with 
logiVm = 19.32 and logAWi < 14.3 which lies 86 kpc 
from a sub-L* galaxy. These four absorption systems, 
with multiphase ratios log [Nm/Novi] > 2.5, 2.9, 5.0, 
and 5.0 respectively, may be shielded from external ion- 
izing flux by adjacent IGM clouds and have unshocked 
gas. 

Mctallicity is one of the great unknowns in interpret- 
ing the baryon content of the WHIM absorbers. In pre- 
vious studies, (O/H) has been assumed to be constant 
at 10% solar. However, one might expect that metallic- 
ity variations could play a role in the multiphase ratio: 
outflows of material from galaxies should be denser and 
more enriched than primordial IGM clouds. This would 
contribute a negative slope in the multiphase ratio plot 
(Figure lb). Instead, the denser (high-Ani) absorbers 
show a relatively lower ratio of O VI/H I, analogous 
to lower metallicity. While metallicity variations almost 
certainly exist in the IGM, they are not the main cause 
of the positive slope in the multiphase ratio plot. Some 
of this effect could arise from changes in the mean O VI 
ionization fraction, typically chosen to be the maximum 
value, fovi ~ 0.2 at T may = 10 5 45 K in collisiq nal ion- 
ization equilibrium ({Sutherland &: Dopital 119931) . This 
fraction is expected to vary, depending on the range of 
shock velocities that prod uce the O VI l)Heckman et, alJ 
120(51 IRaian fc ShulllSool . 

4. THE HOT BARYON CONTENT OF THE UNIVERSE 

Using O VI as a tracer of the 10 5-6 K portion of the 
WHIM, we can employ the number of absorbers per unit 
rcdshift, d/Vovi/ dz, to determine f^wHiM, the fraction of 
the critical density contributed by this WHIM gas. Our 
detection statistics are shown in Figure 2 for O VI. 
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Fig. 2. — Column density histogram for the observed OVI ab- 
sorbers. Approximate equivalent widths for OVI A1032 are shown 
along the top axis. The decreasing number of absorbers at low 
column densities suggests that our survey is incomplete for weaker 
absorption lines. 



The total redshift pathlength of our survey is a func- 
tion of equivalent width, and our survey is more com- 
plete for strong absorbers than for weak absorbers. The 
equivalent- width sensitivity, W m ;„(A), is a function of 
spectrograph resolution, R = A/AA, and the signal-to- 
noise ratio (S/N) of the data per resolution element. 
Thus, for a 4cr detection limit, we define W m i n (A) = 
4{\/R)/{S/N). We assume R = 20,000 and calculate 
W m in profiles for each dataset based on S/N measured ev- 
ery 10 A. Strong instrumental features, ISM absorption, 
and airglow lines are masked out (S/N set to zero). From 
this procedure, we calculate A r m i n (A) for both O VI tran- 
sitions, using curves of growth with 6ovi = 25 km s _1 . 
By moving an absorber doublet along the profile as a 
function of z a b s , we generate a profile N m - m (z). By 
adding up the total path length for each absorber at each 
N m i n (z), we determine the relationship between iVovi 
and Az as shown in Figure 3. With a total high-iVovi 
pathlength Az — 2.21, we are at least 80% complete in 
O VI detection down to logiVovi = 13.4 (W\ = 30 mA 
in the 1032 A line). However, Az and the completeness 
fall off rapidly for weaker absorbers. Dividing Figure 2 
by Figure 3, we obtain the profile of dMovi/dz as a func- 
tion of Novi (Figure 4). To first order, this procedure 
should correct for incompleteness in our O VI survey. 

Our incompleteness-corrected value of the absorber fre- 
quency, dMavi/dz, compares well with previous, un- 
corrected values (Table 1). For O VI absorbers with 
W\ > 50 mA (in A1032), we find (Wovi/dz = 9 ± 2, 
somewhat lower than the value, dMovi/dz = 14*g, found 
bv I Savage et alJ l)2002|) using six O VI absorbers toward 
PG 0953+415. For the weaker O VI absorbers, we find 
dMovi/dz = 17 + 3 (W A > 30 mA) and dNovi/dz = 
19 + 3 (W A > 15 mA) from our total sample of 40_ 6 
measured O VI absorbers. Uncertainties are based on 
single-sided la confidence limits in Pois s on st atistics 
l!Gehrelslll986T) . ITripp. Savage, fc Jenkins! lEfflM found 
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Fig. 3. — Cumulative FUSE redshift pathlength, Az, as function 
of 4<t absorption-line detection limits in O VI column density and 
A1032 equivalent width. Strong lines can be detected in poor- 
quality data. Weak lines can only be detected in the highest-S/N 
observations and thus contribute less total path length. We use 
the Az profile to correct for incompleteness in our O VI survey. 



dMovi/dz > 17 at 90% confidence for Wa > 30 mA, 
based on four abs orbers toward H18 21+643. 

A recent study (|TriDn et al]l2004l) using STIS/E140M 
data finds dAfovi/dz = 23 ± 4 for W A > 30 mA based 
on 44 O VI absorbers. This sample was taken at slightly 
higher redshift (0.12 < z a b s < 0.57) than our sample 
(z a bs ^ 0.15) and shows a higher value of dMavi/dz. 
This is not predicted by simulations, which predict an in- 
creasing WHIM fraction at recent epochs. However, the 
difference in redshift between the two samples is small, 
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Fig. 4. — Number of O VI absorbers per column density bin per 
unit redshift after our first-order completeness correction. Even 
after this correction, the number of O VI absorbers appears to fall 
off rapidly below log./Vovi< 13.4. However, the statistics are poor 
with only five absorbers in the lowest two bins. 
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Fig. 5. — Cumulative eq uivalent width dist ribution for O VI ab- 
sorbers after Figure 10 of Che n et alj 120031) . Our data points 
(filled diamonds) at (O VI A1032) equivalent width detection 
threshold s of 10, 15, 30 , 50, an d 100 mA are slightly lower than val- 
ues fr om ISavage et alj ^2002|) (sq uare), I'lripp. Savage. fe Jenkins! 
jgOOCh (circle), an d ITripp et alj 120041) (triangle). Models from 
IChen et alj 120031) show the simulated distribution at a range of 
metallicities: Z = O.IZq (solid curve), Z = Zq (dotted), and 
an overdensity-dependent metallicity Z = Z(5+f) (dashed) . The 
Z = 0.5 Zq model used by Fang. Brvan. & Canizarcs (2002) (dot- 
dashed) represents a more simplified physical model. 



and the discrepancy in dNovi/dz may be a result of cos- 
mic variance. 

Recent cosmological simulations of the X-ray forest 
predict a distribution of O VI absorbers and provide 
a convenient way to check curr ent simulations with ob- 
served data. IChen et al.l l)2003t) mo del the X-ray fores t 
assuming a ACDM model similar to iDave et alJ {2001), 
both collisional ionization and photoionization from a 
UV background, radiative cooling, and a range of dif- 
ferent metallicities. The cumulative distribution of O VI 
absorbers (dA/ovi / 'cfln[l + z]) down to a minimum equiv- 
alent width is drawn from these models and provides a 
convenient comparison with observed statistics from our 
work and previous surveys (see Figure 5). 

The simulation with Z = 0.1 Zq shows a reasonably 
good match to the observed cumulative distribution of 
O VI absorbers in the local universe. Weaker absorbers 
{W\ < 30 mA) are overpredicted compared with observa- 
tions, but this may be a matter of small-number statis- 
tics or differing detection thresholds between observa- 
tions and simulations. A simulation in which metallicity 
is a function of overdensity fits the data sli ghtly better 
than t he fixed-metallicity model. However. IChen et alJ 
(2003) caution that there is a substantial scatter in the 
simulations, so that the lower-metallicity curves are es- 
sentially indistinguishable . The Z = 0.5 Z Q model of 
iFang. Brvan. & Canizare'sl l)2002fl also fits the observed 
data reasonably well, but it is based on a less physical 
simulation with no photoionization and no radiative cool- 
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TABLE 1 
IGM O VI Absorber Statistics 



and 



n H 



Criteria 


A<ab S 


d/V /dz 


W > 10 mA 


40 


21 ±3 


W > 15 mA 


38 


19 ±3 


W > 30 mA 


35 


17 ±3 


W > 50 mA 


19 


9 ± 2 


W > 100 mA 


6 


3±? 


W > 30 mA 


Id 


23 ±4 


W > 30 mA 


4 


> 17 


W > 50 mA 


6 


14t 9 6 



WHIM 



reference 



" (1.90 x 10 

;20 W-° Az 



- 7 cm- 3 )(l + z) 3 



(2) 



0.0022 ± 0.0003 
0.0021 ± 0.0003 
0.0021 ± 0.0004 
0.0016 ± 0.0005 
0.0008 ± 0.0004 

0.0027 
>0.006 b 
>0.003 



this work 
this work 
this work 
this work 
this work 



a All f^wHlM values have been converted to a consistent set of as- 
sumptions: /ovi = 0.2, Z = 0.1 Z®, Hp = 70/170 km s^ 1 Mpc" 1 , 
and (O/H) = 4.9 X 10" 4 lAllende Pr igg_gl^ ^0mT) . 

b All four OVI absorbers in Tripp. Savage, fc Jenkins! l l200fj) lie 
along the unusually rich H1821+643 sight line. 



The above relations allow us to relate TVhi to the physical 
gas density, n# , needed for the hydrogen photoionization 
correction. Equation (1) is derived by fitting the multi- 
phase correlation in Fig. lb over the approximate range 
10 < <5 H < 300. The scaling constant C = 2.5 ± 0.2 is 
Tripp et al. (2004) the mean multiphase ratio at log Am = 14, and the best- 
Tripp. Savage, fc ^eiMaUgOfX^^ i s a = 0.9 ± 0.1. Equation (2) comes from 
Savage et ak (2 QQ2) CO smological simulations l|Dave et all 11999(1 and relates 
the hydrogen overdensity, Sh, to the H I column density, 
N H i = (10 14 cm- 2 )iVi4. 

From these relations, we can derive a statistical value 
of the (O /H) metallicity from the formula, 



ing. The solar metallicity model of the IGM can clearly 
be ruled out by our observations. Simulations at lower 
metallicities (Z = 0.01 Zq) would be helpful, given the 
wealth of new observational results. 

From dN /dz, we can calcula te the contribution to fib 
from WHIM discussed in lSavage et all l|2~002'l . We 

assume a Hubble constant Ho = 70hro km s" 1 Mpc -1 , 
rather than their value /175, and we make standard as- 
sumptions regarding O VI ionization fraction and (O/H) 
meta llicity <|Tripp. Savage, fe Jenkinsl2000tlSavage et al.l 
2002). We adopt an ionization fraction, /ovi = 0.2 char- 
acteristic of its maximum value in collisional ionization 
equilibrium (CIE) and assume an O/H abundance 10% 
of the solar value. Of course, CIE is almost certainly 
not a valid assumption for the low-density IGM. As the 
infalling gas is shock-heated during structure formation, 
the ionization states O VI, O VII, and O VIII undergo 
transient spikes in abundance, followed by cooling and 
recombination. Recent time-dependent models of the 
ionization, recombin ation, and cooling of shock-heated, 
low-density WHIM l|Raian fc Shulll feOO.^ find that the 
mean, time-averaged O VI ion fraction is (/ovi) =8-35%, 
over a range of initial temperatures 5.4 < logT < 6.2. 
These fractions are compatible (with a factor of two 
spread) with the fiducial value, /ovi = 0.2, based on 
the maximum fraction of O VI at T max = 10 5 45 K, in 
CIE llSutherland fc Dor)italll99l . 

As shown above (Fig. 5), a metallicity Z = 0.1 Z & is 
reasonably consistent with the observed equivalent width 
distribution of O VI absorbers at z ~ 0. However, using 
our sample of O VI and H I absorbers, we can make a 
more direct estimate of the (O/H) metallicity. To do 
so, we must estimate the amount of hydrogen associated 
with each component of the the multiphase system of 
photoionized H I and collisionally ionized O VI. Since 
the H I and O VI absorbers are associated kinematically, 
we assume they share the same metallicity. We then use 
the empirical relations between iVm and overdensity, Sh, 
and between the multiphase ratio (Nm/Noyi) and A/hi- 



(No/Nh) = (N ovl /N m ) x 



./'] 



HI 



OVI 



(3) 



rr.-0.726p-! 
1 A 1 -1? 



We employ the multiphase ratio, Nm/ Novi, together 
with appropriate ionization correction factors, /ovi an d 
/hi- The H I fraction is derived from photoionization 
equilibrium in the low- z IGM 

/hi = ^ = (1.80x10-) ( 1+^(| 

for gas with n e = 1.16n# at temperature (10 K)T4, 
photoionized at rate Th = (10~ 13 s _1 )r_i3. Since 
the mean-free path of a Lyman continuum photon is 
very large in the low-overdensity IGM, we use the 
case-A recombination rate coefficient ci^P — (4.09 x 
10 -13 cm 3 s" 1 )^ -0 ' 726 . Combining the two empirical re- 
lations with the relation of photoionization equilibrium, 
we find that the mean oxygen metallicity of the O VI 
absorbers at (z) = 0.06 is 



Z = {QmZ Q )N^- 2 T^ J2& Tl 



/ovi 
0.2 



(5) 



Given the uncertainties in these empirical relations, it is 
remarkable that this formula arrives at an oxygen abun- 
dance near the fiducial value of 10% solar. In fact, our 
estimated value, Zq = 0.09Zq, is probably accurate to 
only a factor of 2. 
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We now derive the fractional contribution to closure 
density of WHIM baryons associated with hot O VI, 



9. 



WHIM : 



c Pc J (0/H) Q Z /, 



ovi 



N n 



W mln 



(1.85 x 10" 



(Novi)dN OVI 



(6) 



We perform the above sum using the value, {dhf / dz)i, for 
each column-density bin in Figure 4, with {Ncyyiji cho- 
sen as the mean column density (cm -2 ) in the bin. We 
find ftwHiM = (0.0022 ± 0.0003)^ . Values of ^whim 
for other equivalent width thresholds are listed in Ta- 
ble 1 along with the results of other studies converted to 
uniform values of Hq, /ovi, an d (O/H)©. 

It should be stressed that these estimates of SIwhim are 
based on an assumed uniform 10% solar (O/H) metallic- 
ity and that further uncertainty arises in our assumed 
O VI ionization fract ion. Changes in fn yj within the 
range of mean values jjlaian & Shull 2005) could change 
our result by a factor of two. Based on the standard 
assumptions, our result is s lightly lower than the result 
published by iSavage et alJ l)2002lh fiwfflM > 0.002/if,; 1 . 
iSavaee et alJ l)2002j) used an older solar oxygen abun- 
dan ce, (Q/H) m = 7.41 x 10" 4 which is 50% larger than 
the lAllende Prieto et all l|200j) value. After correcting 
for this difference, their result for fiwfflM is larger than 
ours. iTripp et all l)2004[) find ^whim = 0.0027 based on 
44 absorbers at z > 0.12. 

Our estimate corresponds to f2wHiivt/^& = 0.048 ± 
0.007, so that WHIM gas in the range 10 5 ~ 6 K makes 
up roughly 5% of the baryonic mass in the local uni- 
verse. The contribution could exceed 10% if we ac- 
count for the likely spread in (O/H) metallicities and 
O VI ionization fractions. Simulations predict that the 
current universe is composed of ~ 30% WHIM gas, so 
our value of ~ 3 — 10% falls short of this mark. How- 
ever, O VI is only useful as a proxy for gas within the 
lower portion of the WHIM temperature range. Account- 
ing for the hotter 10 6-7 K gas using similar method- 
ology will require high-sensitivity, high resolution X- 
ray observations of O VII and O VIII with spectro- 
graphs ab oard future missions such as Constellation-X 
orXEUS l|Fang. Brvan. fc Canizaresl l200l IChen et all 
|2003(). Limited observatio nal work has been done in this 
area ijNicastro et al.ll2005|) , but the current generation of 
X-ray telescopes are not ideal for this kind of investiga- 
tion. 

The distribution of H I absorbers with column 
density is often expressed as a power law w i th in- 
dex 0: dMIdN oc JV-f llWevmann et all 119981 
IPenton. Stocke. fc S hull 2000, 2004). 



For the first time, 



this analysis can be applied to WHIM species. We find 
that the differential number of O VI absorbers with col- 
umn density is a power law, dMovi/dNovi oc N^yf^ ' 1 
for absorbers with logAfovi > 13.4. This is some- 
what ste eper than the corresponding H I distribution, 
~ 1.6 l|Penton. Stocke. fc ShuHl2004) and means that 
low-column absorbers contribute comparable amounts to 
the O VI baryon census as the rare high-column systems. 
Note that in eq. [6], ^whim scales as N~° n 2 for = 2.2. 
Even though we correct for incompleteness in the weaker 
absorbers, the statistics are still poor at log iVovi < 13.4. 
It is unclear whether the turnover in dAfovi/dz at lower 
columns is real or a statistical fluctuation from small 
numbers of weak absorbers. Nevertheless, the weak O VI 
absorbers appear to make significant contributions to the 
baryon mass density. The steep power law also reinforces 
our conclusions about the nature of the multiphase IGM 
with a core-halo structure. 

We have demonstrated that the contributions from 
weak O VI absorbers cannot be neglected in an accurate 
WHIM baryon census. The rare strong O VI absorbers 
and numerous weak absorbers contribute nearly equally 
to fiwHiM- Further FUSE observations at high-S/N will 
allow us to probe weak O VI absorbers and refine the 
statistics at the low-column end of the absorber distri- 
bution. An analysis of O VI absorbers at higher redshifts 
[z > 0.5) would allow us to track changes in the WHIM 
density, confirming the decrease in the amount of shocked 
gas at higher redshifts predicted by cosmological s imula- 
tions ijCen & Ostrikerll 19991 iDave et alJll999l 120011) . 

Only by enlarging our sample of multiphase (H I, O VI) 
absorbers beyond the 40 discussed here will we be able 
to refine our statistical estimate for the metallicity of the 
WHIM. With a much larger sample of O VI absorbers, we 
could use maximum-likelihood techniques to search for 
the expected tr ends: (1) decreasing met allicity at lower 
overdensity Sh l|Cnedin fc - Q j s^rikeJl i 997) : (2) decreasing 
ilovi at higher redshift (Da ve et al.ll999|) : and (3) trends 
of Nm/Novi with other metal indicators such as C III, 
C IV, or Si IV. 



We are grateful to Steve Penton, John Stocke, Jessica 
Rosenberg, Todd Tripp, Bill Blair, Blair Savage, Ken 
Sembach, and Jason Tumlinson for useful discussions re- 
garding this project. This work contains data obtained 
for the Guaranteed Time Team by the NASA-CNES- 
CSA FUSE mission operated by the Johns Hopkins Uni- 
versity, as well as data from the Hubble Space Telescope. 
Financial support to the University of Colorado has 
been provided by NASA/ FUSE contract NAS5-32985 
and grant NAG5-13004, by our HST Lya survey pro- 
gram 6593, and by theoretical grants from NASA/LTSA 
(NAG5-7262) and NSF (AST02-06042). 



REFERENCES 



Allonde Prieto. C, Lambert. D. L.. & Asplund, M. 2001. Ap.J. 55(5. 
L63 

Cen, R., & Ostriker, J. P. 1999, ApJ, 519, L109 
Chen, H-W., Lanzetta, K. M., Webb, J. K., & Barcons, X., 1998, 
ApJ, 498, 77 

Chen, X., Weinberg, D. H., Katz, N., & Dave, R. 2003, ApJ, 594, 
42 

Collins, J., Shull, J. M., & Giroux, M. 2004, ApJ, 605, 216 



Danforth, C. W., Shull, J. M., & Rosenberg, J. L., 2005, ApJ, in 

preparation (Paper II) 
Dave, R., et al. 1999, ApJ, 511, 521 
Dave, R., et al. 2001, ApJ, 552, 473 

Fang, T., Bryan, G. L., & Canizares, C. R. 2002, ApJ, 564, 604 
Gehrels, N. 1986, ApJ, 303, 336 

Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, 581 



7 



Grcvcsse, N., Noels, A., & Sauval, A. J. 1996, in ASP Conf. Ser. 

99, Cosmic Abundances, ed. S. S. Holt & G. Sonneborn (San 

Francisco: ASP), 117 
Heckman, T. M., Norman, C. A., Strickland, D. K., & Sembach, 

K. R. 2002, ApJ, 577, 691 
Jenkins, E. B., Bowcn, D. V., Tripp, T. M., Sembach, K. R., 

Leighly, K. M., Halpern, J. P., & Lauroesch, J. T. 2003, AJ, 

125, 2824 

Lanzetta, K. M., et al. 1995, ApJ, 442, 538 
Nicastro, F., ct al. 2005, Nature, in press 

Pcnton, S. V., Shull, J. M., & Stocke, J. T. 2000, ApJ, 544, 150 
Pcnton, S. V., Stocke, J. T., & Shull, J. M. 2000, ApJS, 130, 121 
Pcnton, S. V., Stocke, J. T., & Shull, J. M. 2003, ApJ, 565, 720 
Penton, S. V., Stocke, J. T., & Shull, J. M. 2004, ApJS, 152, 29 
Rajan, N. & Shull, J. M., 2005, ApJ, in preparation 
Savage, B. D., Sembach, K. R., Tripp, T. M., & Richter, P. 2002, 
ApJ, 564, 631 



Savage, B. D., & Sembach, K. R. 1991, ApJ, 379, 245 

Sembach, K. R., ct al. 2003, ApJS, 146, 165 

Shull, J. M., Giroux, M., et al. 2000, ApJ, 538, L13 

Shull, J. M. 2003, in The IGM/Galaxy Connection: The 

Distribution of Baryons at z = 0, ed. J. L. Rosenberg & M. 

E. Putman, (Dordrecht: Kluwcr), 1 
Stocke, J. T., et al. 2004, ApJ, 609, 94 
Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253 
Tripp, T. M., Savage, B. D., & Jenkins, E. B. 2000, ApJ, 534, LI 
Tripp, T. M., et al. 2004, in Astrophysics in the Far Ultraviolet: 

Five Years of Discovery with FUSE, Victoria, BC. 
Tripp, T. M., et al. 2005, ApJ, in press 
Weymann, R., et al. 1998, ApJ, 506, 1 



